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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

TECHNICAL MEMORANDUM NO . 1233 

CONTRIBUTION TO THE PROBLEM OF FLOW AT HIGH SPEED* 
By C. Schmieden and K. H. KavalM 

OUTLINE 


Part I. A Few General Remarks Covering the Prandtl-Busemann Method 

Part II . Effect of Compressibility in Axially Symmetrical Flow around 
an E lli psoid 


PRELIMINARY REMARES 


The authors regret that due to lack of time the following investiga- 
tions could not be carried out to a more finished form. Especially in 
the first part it was intended to include a few further applications and 
to use them in the general considerations of this part. In spite of the 
fact that the intentions of the authors could not be realized, the 
authors felt that it would serve the alms of the competition to present 
part I in its present fragmentary form. 


*Beitrage zum Umstromungsproblem. bei hohen Geschwindigkeiten. " 
Lilienthal-Gesellsohaft fur Luftfahrtforschung Bericht S 13 /l* Teil, 
19l2, pp. 10-68. (Figures referred to in Part I are found immediately 
after the appendix to that part.) 
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PART I. A PEW GENERAL REMARKS 
COVERING THE PRANDTIr- 
BUSEMANN METHOD 

INTRODUCTION 


For the solution of the problem of flow in the subsonic range, 
two approximation processes are available, the Janzen— Rayleigh method 
which, proceeding from the potential or the stream function, represents 
the velocities in form of a progressive series of powers of the Mach 
number, and the Prandtl—Busemann method which, based upon the fundamental 
concept of the Prandtl rule, determines the velocities by an expansion 
in series according to a geometric parameter characterizing the body 
(references 1 and 2). Both methods supplement each other most opportunely 
insofar as the first is suited particularly for thick bodies, hence for 
relatively small critical Mach numbers, while the second method works 
best for slender body contours whose critical Mach numbers are close to 
unity. Little is known so far about the limit of convergence of the 
employed series expansion of both methods; no proof has as yet been 
given concerning the ceasing of convergence upon the reaching of local 
sonic velocity; however, it appears plausible from the construction of 
the, at times very small, number of the explicitly computed terms of 
the series. The first section of the present report utilizes the long- 
known relation between velocity and stream density (pV) in order to 
prove that with the use of the stream function the differential equations, 
which serve as basis for the practical calculation in both methods, lose 
their meaning if the local sonic velocity is exceeded anywhere in the 
field of the flow. The fact that the solutions of these differential 
equations obtained by iteration then become useless also, is not sur- 
prising. In section 2 it is shown that the expression employed by 
Hantzche and Wendt for the stream function is not general enough and 
therefore is likely to fail under certain circumstances , as for instance, 
for the flow around an ellipse. The latter is discussed in detail and 
on the basis of the results obtained the suppositions are given and made 
plausible concerning the fact that the series expansion of the velocities 
computed by the Janzen-Rayleigh, or the Prandtl—Busemann method can be 
made to agree by formal transformation. Finally, in the concluding 
section, the arguments of the first two sections are evaluated for the 
problem of convergence of the series expansion built upon the potential, 
with special reference to a report by Gtfrtler (reference 6). 
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1 . If one desires to treat a compressible flow with the differential 
equation 

^(a 2 - u 2 ) + tyy(a 2 - v 2 ) - St^uv = 0 (l) 


of the corresponding stream function In the two-dimensional case, the 
quantities a, u, v must be represented as functions of the 
derivatives of so as to obtain a differential equation containing 
only and its derivatives. In the analogous case of the potential 
this transformation is very simple, since owing to 


u a 


V = 


4 > 


7 


one immediately obtains 


^(a 2 - 3» x 2 ) + <Djry(a 2 - <f y 2 ) - 24^$^ = 0 

> 



(la) 


whereby this equation is obviously valid for all physically possible 
velocities. But the use of the stream function in this instance involves 
a characteristic difficulty which becomes most readily apparent through 
elimination of the sonic velocity from equation (l). If n is the direc- 
tion of the normal in a point of a stream line, the relation 
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>1 


exists ‘between stream function ^ , velocity Y , and stream 
density pY, while for adiabatic flow 


£_ = 
Po 


1 


Lll ;'V 


,'H * 1 

Vj 


V 


1 

K-i 


( 3 ) 


is valid* 

From ( 3 ) it is seen that p/pg decreases monotonic ally, with 
increasing velocity) substituting (3) into (2) it is apparent 
that, while D increases at first with increasing velocity, that 
on passing through the critical velocity - V = a* - it reaches a 
maximum and in the supersonic range for Y approaching V max it 

decreases to zero again- The approximate course of D against V 
is shown in figure 1. To effect tlje required el imina tion of a, 

V is used as function of D or that is, the inverse of 

on 

the above function D(V), which will be, in order to be able to 
make practical calculations, the inverse of the form of a power 
series 


V * £(D) 




The problem is thus the following: 

Given: a function D(v) 

Required: the inverse of this function in form of a power 
series (if), such that V = 0 for D = 0- 

This problem, as known, is solved by the Lagrange inversion formula. 
This Lagrange series converges in a circular disk with the radius 
I^axJ in this case equals the value of D, for which 

= 0. This result, which is also clear from inspection of the 
dV 

sketch, naturally holds also for the case where the inverse function 
is not developed at D = 0 but rather at' some other point between 
0 and L&ax> as is &one in the Prandtl-Busemann method- The 
radius of convergence is then correspondingly reduced such that 
the limit of convergence coincides with sonic boundary. The same 

S\l/ 

holds for the inverse functions u and v as functions of 

ox 
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and vhlch are necessary for the transformation of (l), 

ay 

although in this case the relatione are not as clear as 
for sonic velocity* Xf nra substitutes such power developments 
in (l) the validity of the thus obtained differential equation 
for ^ ceases on reaching sonic velocity. Consequently, one 
should not expect the iterative solutions of the differential 
equations to have significance once the sonic boundary is 
exceeded in them. The sane applies to the RaylaSgh method vhen it 
uses the stream function, as well as to the Hantzeche -Wendt 
method, (developed on 'the suggestion of Busemann) which improves 
the Prandtl rule by successive approximation. Both methods use 
the expansion in. powers of the velocity components, hence must 
fail on reaching the sonic boundary. 


It should be interesting to study this failure for an actual 
single case. Consider the flow around the ellipse at zero 
incidence for which the Busemann correction is computed in 
appendix I, equation (20). Even if the stream function computed there 
did retain its validity when a supersonic zone, no matter how- 
small, has developed about the end point of the small axis, 

the curve obtained by plotting ^ against the circumference of 

the ellipse would have to show a course as sketched in figure 2, 
because of the relation between V and 2x, sketched in figure 1. 

Sty 

The quantity ^ would have to attain its maximum at the beginning 


of the supersonic zone and then decrease monotonously to the end 

point of the small axis* Actually ^ shows an entirely different 

course (fig. 3)* It increases monotonically with the arc length 
for a supercritical Mach number and numerically remains always 


smaller than 


But to such a 


max 


3 ^ 


distribution there 


belongs, according to figure 2, for reasons of continuity, either 
a pure subsonic or a pure supersonic flow, both of which are 
impossible at a supercritical Mach number and in the presence 
of stagnation points. 


In figure 
for an ellipse 


k the streamr-density distribution 



distribution^ 


of axis ratio l/lO at three different Mach numbers 
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is plotted against the conventional parameter angle that is 
only the range of str eamr-dens i ty Talues of interest here. It is 
seen that the curve for the subcritical Mach number 0.80 runs below 
the curve for 0.85, in the vicinity of the end point of the small axis, 
hence, that the stream density increases as yet with rising Mach number. 
M = O.85 is almost exactly the critical Mach number, consequently, 
the solution must fail at still larger Mach numbers. As a matter of 
fact, there results for M = 0.90 the above sketched behavior of 

the ^ distribution; the solution is impractical, 
da 

Figure 5 serves to illustrate the good convergence of the 
Busemann method in the range of its validity. For the same ellipse 
as in figure 4 the stream density at the end point of the minor axis 
is again plotted as function of the Mach number in the range of 
interest; the upper curve holds for the first step of the method, 
the Prandtl rule, the lower curve for the second approximation. The 
intersection of these curves with the straight line of the sonic 
boundary gives the critical Mach number in the corresponding 
approximation (hence M cr = 0.74 and O.78). 


In connection with the foregoing arguments, however, another 
value is of greater Interest, -which can also be taken from these 
curves. The maximum stream density in first approximation has the 


_ PV f pV 4- 

value 0.598 -■■■■■ in the second approximation the value O.583 
p o a o 


p o a o 


whereas the true value naturally coincides with the sonic boundary. 


pV + 

hence must be O.578 ®ius, since the rigorous solution should 

p o a o 

yield a stream-density distribution which would have to approach the 
sonic boundary from, below, at a critical Mach number of around 0.6l, the 
plot Indicates that the difference between the first approximation and 
the rigorous solution is largely canceled out by the second step of the 
method, so that the method converges quite well. 


Moreover, it is pointed out that, at variance with the customary 
terminology in the discussion of figure 4, the Mach number used as 
critical Mach number is that corresponding to the maximum of the stream 
density; for up to that Mach number the solution obtained by this 
approximation has a meaning. 



NACA TM No. 1233 


7 


2. Hi© Hantzsche-Wendt calculations proceed from the following 
assumption. If X is a parameter that characterizes the departure of 
the profile form from a straight section profile (for example, thick— 
ness, camber, angle of attack) at zero incidence the stream function 
can be represented by a power series in powers of X, the coefficients 
of which are functions of x and y: 

■yjr p 

jf 3 y + e(*ir) * + i»(*y) k + . . . (1) 


On the basis of this assumption it is possible to set up, by 
means of comparison of coefficients, a recursive system of differential 
equations for the functions of the individual powers of X which then 
can be integrated as in the Janzen-Rayleigh method. 

But while in the Rayleigh method the expansion of the stream 
function in powers of 


♦ 






( 2 ) 


is always possible, the expansion of f in the form (l) is by no 
means guaranteed a priori. On the contrary it will be shown that it 
is not possible at least in the two cases cited here, the ellipse and 
the ellipsoid of revolution both at zero incidence^ but that sooner or 
later terms of the form In X must appear in these cases. If 
that is the case it ceases to be possible to set up the recursive system 
of equations as was done by Hantzsche— Wendt. Hie best procedure is 
then as follows: The stream function is put in the form 

u = V^ x i y ^ + ^i( x iy) + * 2 ( x i y ^ + • • • (3) 

Tjr 0 is the stream function in the undisturbed free stream, to + tp 
represents the Prandtl approximation; the functions to with u > 2 

which are assumed to satisfy the boundary condition t y = 0 c>n 
the contour, are 'regular in the outside d omain and at infinity 
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their derivatives vanish to a sufficient degree. It is further 
assumed that, apart from sufficiently small zones, especially 
the iamediate vicinity of stagnation points, the function tyy 
and its derivatives are small compared to all and their 
derivatives with subscripts K < u, The question whether this 
assumption is satisfied can be answered only after calculation 
of 

The fact that the concent of the order of magnitude is a 
little vague, especially since the vicinity of the stagnation points must 
be excluded, lies in the nature of the problem and seems therefore 
unavoidable; besides no difficulties arose in the calculation 
of actual cases. When limited to the second approximation the 
new rule gives the same result in all the cases treated by 
Hantzsche -Wendt, in particular, the expansion parameters 
assumed in (l) are automatically obtained- But if 1/ can no 
longer be represented by a power development of the form (l) , 
the differently constructed terms are covered by our formula. 

The necessity of excluding the stagnation point vicinity from 
the appraisal of the order of magnitude lies in the initial Btep 
of the method, the Prandtl rule, and has been voiced often enough 
as principal objection against this rule when stagnation points 
exist. More accurate numerical checks (reference 3) have shown, 
however, that the Prandtl rule gives a fairly close approximation 
even in the vicinity of the stagnation points, although its 
assumptions are by no means satisfied any longer. It is not 
believed that thiB behavior is due to a lucky accident. However, 
the following supposition (unfortunately without proof) seems 
reasonable : 

If it were possible to compute the series (2) to any high U, 
this series would in all probability converge absolutely, but 
not uniformly and would merely represent a rearrangement of 
the Rayleigh series computed for the same contour to any high 
powers of Mach number. 

In other words, if each tern in (2) is developed in powers 
of the Mach number the formal rearrangement in series of these 
powers gives the Rayleigh series. 

If this supposition is correct, the Prandtl rule would be 
legitimate also for the vicinity of the stagnation point in a 
certain manner and would explain the surprisingly good approxi- 
mation of the rule at these points. Unfortunately there is very 
little material available to test this conjecture* whicli for 
lack of time could be evaluated only in two cases, the ellipse 
and the circular-arc profile with shock-free entry. 
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In both cases the result of the examination was poeitivoj 
however, only the more interesting case of the ellipse is treated 
here* 


The Rayleigh approximation up to the term with inclusive 
is known (reference h ) . The end point of the small axis, that is, 
the maximum speed is 


= — 1 — J g + fl_jL_o£ _ £ ^.l.±J. £ I g 2ji ^±jj 

u rr^\ 2 I 0 2 i . «r 


+ € 


3 1 + 2a _ it 


o3 


arc tan 


1 - cr 2 


l 


i.“ o g 
1 + o 2 


(4) 


The development in powers of e gives, up to the terms with £ 


I-l + e + 

U 2 


€ + e c + € 


Hi' 1 ) 


+ c 4 In c 


(5) 


On the other hand the second approximation computed by (l) or (3) 
gives in both cases, if one develops in terms of Mach number and 
breaks off with the term in M 2 , 

3sg£ 3 i + € + £. ( € + € a) (6) 

so that the terms, which are for e as well as for M are at 
the most of the second degree, will be contained in both formulas. 
The differences in the definition of the Mach number - Rayleigh 
ref ere to the velocity of sound at rest formulas (l) and (3) to 
the free-stream sonic velocity - do not make themselves felt as 
yet in this approximation* 

The third approximation in the appendix of the Busemann 
method based on (3) was carried only far enough to show the 
appearance of a term In e* It was not further carried out 

because of the prohibitive amount of calculations required in order 
to obtain the final formulas from the final formula given there; 
it can be seen that a term - e - £ must appear too, vhich 



10 


NACA TM No . 1233 


likewise implies that agreement of both developments is to be expected 
for the terms In e and e3. 

3 . The following conclusions are drawn with the use of the stream 
function. The Prendtl— Busemann and the Janzen— Rayleigh method give 
identical results, but only when both developments are carried out 
completely. Each development represents a rearrangement of the other. 

Both developments diverge if at any point in the flow domain the local 
velocity of sound is exceeded. 

But in both cases the velocity potential could be used instead of 
the stream function, as is, in fact, done predominantly in the Rayleigh 
method. However, no case of flow past a body has been calculated as 
yet with the Prandtl-Busemann method. Ctae of the likely reasons for 
this might be found in the complicated nature of the boundary condition 
for the potential. Since the reduction of the differential equation to 
coefficients dependent on 4 alone is in this instance easily achieved 
for the entire physically possible speed range (as mentioned already 
in (la)), no absolutely valid reason for a ceasing of convergence of 
the solution on passing through the velocity of sound can be found 
from the differential equation alone, owing to its nonlinearity. But, 
since the values of the velocity computed for identical conditions 
from the stream function and the potential are identical in the domain 
of convergence of the stream function, there is a strong suspicion that 
the potential development itself ceases to converge at the same place. 

In faot, all the calculation made by the Rayleigh method with the potential, 
particularly those by Lamia (reference 5) for circle and sphere, indicate 
that the obtained series cease to converge on reaching the sonic boundary, 
so far as ouch a conclusion can be drawn at all from the few explicitly 
known terms of these series. Consolidating this result with those 
obtained above for the stream function it may be stated with great 
probability that even the fourth method, -namely, the Prandtl-Busemann 
method which uses the potential, has the same limit of convergence a 3 the 
other three. This statement contradicts the result of a report by 
GSrtlor (reference 6 ); therefore, a brief critical review of his report 
will be given consequently. 

GoVtler used his method to compute the flow past a wavy wall where 
supersonic zone of finite extent occurs. Against this result of Gdrtler 
some serious objections may be made. 
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1. The problem of convergence of his development, of which 
the first three terms are explicitly calculated, is practically 
ignored. 


2. A part of the coefficients of hie first example is wrong, 
according to a careful check and with the correct coefficients 
the convergence becomes definitely worse. 

3* The curves of figures l(a) and l(b) ill his report are 
constructed with the wrong coefficients. T/hether the correct 
values of the constants already bring the supersonic zone, small 
by itself, to disappearance could not be checked, unfortunately. 

This is not the place to be polemical concerning GCrtler's 
report. That was not our intention. However, the doubts have to 
be brought to attention, for his report does not give a proof, 
that with the development followed therein a solution could be 
obtained, whose range of validity would extend to a small region 
of the supersonic domain. 
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APPENDIX 

OAJjCULATIC'N OF THE THIRD APPRCtSIMAIlCSN PCR THE ELLIPSE 
In the general differential equation for the stream function 

'J'jjjX* 2 " u 2 ) +^yy(a 2 - v 2 ) - « 0 (l) 

it is necessary, in accordance vith the Prandtl-Buseaenn method of 
the splitting off of the free -stream 1 part U parallel to- the 
r-axis, to set 


U *1 U 4- u^) V » V^ 


( 2 ) 


where the quantities carrying the subscript 1 denote disturbance 
quantities. The following is then rigorously valid (a^ = velocity 
of sound of free stream U) 


a 2 a Oj 2 - (k - 1) U u^ - " (u^ 2 + V^ 2 ) 

u 2 « TJ 2 + 2U Uj_ + u^ 2 j v 2 c v^ 2 
We write £or abbreviation 

M * U » — rst:?*! -* |i 2 M 2 * U 2 - 1 (4) 

°1 \/l - M 2 

Equation (3) substituted in (l) gives 

^[“l 2 - B 2 - O' + 1) O uj. - £-±-i Uj 2 - JS-g-i | 

^yy[^ - O' * 1) 0 Uj, - Lfl V^] 


-S^XyftJ Vj_ + Ul V-J_) = 0 


( 5 ) 
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The abbreviation X * u * i gives for the subsonic son© 

Po * 

in second approximation the relation between the disturbance 
velocities and the disturbance stream function 1 ]/^. (Owing to 

^xx “^lxx et0 ' tlle 8U ks cri P t 1 for is °nii t ' fc © i hereafter 

ae there exists no danger of confusion.) 

r _ 

^% y + B^ .^ y g (k + 1) + U 2 (2 - k )\ + xjTj 2 ‘ 


Y » - S^Ji 2 ^ 


V 


x Y y 


Ui 2 =U 2 U^y 2 
v- 2 = y 2 ^. 2 


(6) 


Dividing (5) by and moving all terms of higher order to 

the right gives 


*„(! -ift ■‘•V *'*'=£ 


( K . , v Hji + Jfi-4=. 3- + S 1 It 

(K + 1) ^g- + 2 ^ + ^2 




yy 


, ^ u u i , K - 1 ^ + — •h-J. Ilf. 

\K - 1 ; 5 — ♦ « o ^ n o 


8 ^ 8 *1 8 


+ 2if ( Lh + 3u!n ( 7) 

“yUj. 2 , 2 U7J 


(6) substituted in (7) gives 


*« <l ’ A * V - ^ 


(k + 1) + (K - 1) U 2 ^y >Ir 


yy 


" 2 *x V 


-u 2 (2 -K) 


+ v. 


+ ijr 2 
x 


yg I 
« 2,^2 ^ 


(K + 1 ) \J/ 2 

y 


-T 




y' 


yy 2r 




si 


U 6 (K + 1) - 2^(K - X) 

ll 


(K + 1) (n 2 - 1) + K - 1 | 

(K - 1) \i/ y 2 [u 6 (K + 1) - 2n 4 (ft - 1) - n 2 (2 - K)j 


-^ 2 2 | ( K " 1) (u 2 " 1) 


l] 

+ K + 1 i > - 






J. 2 V^^y~ (2^“1) 

J 1 


(8) 
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The flow pattern is now affinely distorted in the direction normal 
to the flow by the Prandtl transformation 


i * I; y - J V^yy " 'W 'I'xy = J** 




we further put 


■ X U n v|/* 


and immediately omit the asterisk for after which (8) gives 
by multiplication with n* 


^ c (^ 2 “ 1 ) [(k + 1 ) + (* • D^n 

- 2 'I'g'l'gj + | (n 2 - 1) ^* 2 + 1) \J/ n 2 ^ (K + 1) 

'2 H 2 (k - 1) - (2 - k)] + 'I'g 2 [n 2 ( K + 1) -fijf 

f p w — 

4 nn < (k - 1) tyjj 2 + l) - 2n 2 (K - 1) - (2 - K)J 




2 


[^ 2 (vc - 1) + a]|- 4 ^^ ^ j2p2 - 1^ 


Next we put "'l' i + ^2 + ^3 ♦ ■ • • determine vjj^, \|f 2 , \{/^ ... 

eucceseively from ( 9 ) on the assumption, that \[f is small 'compared 
to (free stream) tyg snail relative to \{/_, small 

relative to 2 and so on. If the flow involves a body with 

stagnation points, this assumption is certainly not fulfilled 
in the vicinity of the stagnation points; in such a case it 
must suffice to postulate the validity of the assumption outside 
of a certain range from the stagnation point and to check the 
Justification for the assumption on hand of the result of the 
calculation. Discounting this difficulty for the present, the 
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conventional iteration method can he used in our case as it was 
used in the Rayleigh method. Thus it obtains for and \J/g 

the following two differential equations 


A^ = OJ A ^ = fa 2 - 1) 

j> + x) •* a +i,+i5g + < K ’ - 6 ’•'ig+jg^ 


(9a) 


Substituting in the remaining terms of the right-hand side 
of (9) gives a part of the right-hand side for the rest of 

the disturbance terms belonging to BTe computed later- 

Because 3 the first part of the disturbance terms 

becomes 



Because A\J/^ » 0, can be regarded as imaginary part of an 

analytical function w(£) £ a (l + in): 


V 


t+ij 


w 


*iT, + i 


t 


iS 


V 


rr 



* *Jgil + 1 *lgg 


On denoting conjugate complex q,H 0 Btities by a bar, all the combina- 
tions of derivatives of \Jz-j_ appearing above can be represented 
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"by I m aginary parts of functions of 


£ and f : 


2 ^ ^ W ” W> + W " W *^ 

2 ^ 9 J (v" v' - w" w 1 ) 

1^) ^4 

" J L w " + 

^ = ’ J 0 ^ w ' ’ *^ 2 ] 
4 ^ig» - j t* < w ' 2 - * , 2 i) 


(ii) 


Substituting (11) in (10) gives 

A 1 ^3 “ J (A v" w * 2 + 2 B v" w' 5' + C w" v' 2 ) (12) 

where 

A » | (n 2 - 1) (a - |J - 7 ), B = | (n 2 - 1) (a + £) 

C » | (u 2 - 1) (a - 3 + 7 ) 

and 

a = + l) n 2 - (k - 1 )"] [n* ( K + l) - 2 (i 2 (k - 1) - ( 2 - K j] , 

3 = u 2 (u 2 - !),(« + 1)' - 2 


7 » b (2 h 2 - 1 ) 
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With the abbreviation 

C 9 ( K + l) (U 2 " 1) 

th© expressions for ©t* 0, y end A, Bj C can b© further consolidated 
a • (o + 2) |jCm- 2 - l) a + ^ “ 3~] 

0 « M 2 a - 2j 7 m k ((? a 2 - 1). 


hence 


S 


(n 2 . 


» - i 

C * - G* 2 
8 


l) 2 a (c + 5) 

1) jj[u 2 - 1) cr 2 + (7 M- 2 
1) {fc 2 - 1) 


a 2 + 5 (h 2 


5) o + 8 (u 2 - l)j >(13) 
1) a + 8 (? a 2 - 3)j 


These terms are now supplemented by those obtained from the products 
of \{/ 1 and after insertion in the first bracket on the 

right-hand side of (9)} they are 


(w 2 - l) 


* (k - l) (t ln * * 2 r, Vi) 

- 2 (v^ + ^v)j 

- (d 2 - i) | ^2n ^1^2, |(k + i) n 2 - (k - i)] 

V _ 

- 2 +2g^ n + *l n [< K + X > “ 8,f '2 S g + < B • 15 *2n, 

- 2 i (14) 
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Now 

^2 * J Q (see also (20)) 
and 

Wg " W£ * Vj; W n = i (wj - wj)j \g n = l (W££ - wff) 

Wgg=W^+ 2 Vf E { + wjj; - - W tt + 2 > £ ; - VjJ 

Hence 

J (WJ + ° J (1 w s " 1 w c } 

+«„ - J <* M n - 1 W S£> J *agg • J < w u ♦ 8 w tr ♦ w f? 

'*'2-1-1 “ J <- w » ♦ 2 W U • *1$ 

and 

2 ^ «v'+v'; 2ii.*v' *5' 

It] A S 

2ty . a v" + j 2 i t. * w" * w" 

1§1 igg 

Substituting (15) and (15a) in (l^) gives 
/• 

A, * **—-1 J - [(k + 1) M 2 - (K - 1)J (Wj - W £) (v" - V") 

- 2 (W£ + W£) (w" + $") + j(K 4- 1) n 2 - (k - l)|(W$£ + Wj£) 

<v' + W r ) - 2 (W ££ - W£g (w' - ff') - 
+ 2 j (k + l) M 2 + (k - 1)J W££ (w 1 + ^’)j ( 16 ) 
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where it should he noted that (see also further "below) 


kv tr a 


~ o v" w‘ + (a + Ij-) w” tf'J 


The right-hand side of (l 6) can be further simplified by the 
previously used abbreviation a 

(k + 1) (h 2 - 1) ** aj (k + l) y. 2 - (k - 1) - o + 2 

, v 2 

<k * 1 ) P +K-l*a + 2 x 


and with it the following representation is obtained for ^ 

/• 

Ag ^3 “ J jcr (W£ w” + v" + W^v* +W^T 7 ') 

- (cr + 4) (Wj V" + w" - ~ w') 

+ tt — I_i (<j + 2k) (w 1 + v') I cr w' w" + (a + 4) w" v'J * 

Equations (12) and (l8) together give the disturbance equation for 
the third approximation ^ in the sense of our iteration method- 

If z Is the coordinate of the image plane which contains a 
circle, £ the coordinate of the plane of flow, we get with R 
as radius of the circle 

_2 ~2 .2 2 v ,2 _2 I 


o 9 2 P 2 

£ = z + — } w = - — ~ < lj e = & — i— =■ 

Z Z R 2 R 2 + a 2 


« £ 


2 _ »2 


dw _ _ R' 


2 . c.2 


dw _ w i = _ Rj 


42 2 2 ' 42 *2 ' 4 £ -2 - *2 
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The transformation to ii&age -plane circle coordinates gives for ^2 
and ¥ 



To reduce the right-hand side of (12) and (l 8) to image-plane circle 
coordinates it ia borne in mind that 
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Hence, after introducing from (12) and (l3) 


aim = a A. (kl) 41 +B f (w' 2 ) 5ut („■) (m'f f 
bz oz dz y 3 J dz <3- z dz o.z [d£ J dE 

U 2 - i f ( a fd_ L' w, ds \5I + d_ (V W,- 5 ^ d$| 

4 H L dz v lz d £ ) dz dz ' d £ / “> 

- (c + 4) [j= W w lB + d_ z v ) Wl _ - d v lz -~'g 


dz 


- d_ [ w,— — \ + 2 (c + 2 k) < a t fx~. lN 

dz V lz d£ J dzj [dz \ 3 / 

+ du.(*!!^ffil + (a + 4)1 d-(3t£\2S + d- ( v >) ^ v 
dz \ 2 / dz ' [dz V. 2 J dz c - z ^- z _ 


It, is essential to note that expressions such as v* W lz d| 

are functions of z and "Z, hence must "be Bplit up before 
(21) can he integrated in the usual manner. 

The function can he written (see also (20) ) 

Wi = Oi (z) I + $ 2 (z) w + (z) 



with 


Hence 


0 ]L (z) - | v ,2 ‘, ® 2 ( z ) * ( c + ^ w< 




(*) 


- & (o’ + *0 W p - 2 -~-k 2 


V 


E c 


W. 


lz 


Y 


£ + V **® 3 ' (z) } W K = 0 X £i * 2 g 
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and 


w n | - 4 i 


+ $ 2 tf'j 



$2 



f 



Correspondingly the explicit expressions for the derivatives ur^d 
in equation ( 21 ) * 


h (*' »u $)& - i *<"'>! S ? £ + '•*’> * « £ 


dz ( W ' w lz df) to * to ("'*1' tf to 

+ fc( w ’ v s*)* g ♦ £ (*• v $) £ 

^ & <*’> *®a‘ *S^*3» 


dzV dy dz dz 2 d y dz + dsH dJdi 


7 ? ( W ir Jf 1 T 2 - (* + 4 ) & 

dz y 12 d^y dz -dz dz d£ 

w,- dt'.* dsL'il + $ dwl & * o &£? ££ 

z dz 1 dz dz 2 dz dz 0 da dz 


+ S_±Ji ds /2 ^ . d_ / &A ^ 

2 dz dz dz y dz dg ) dz 


b \r 
d £ 


(22) 



NACA !WNo. 1233 


23 


The following integrals must he calculated: 

^&£azi is-, /[if'atj Um'iiy/ imu 

dz d£ / dz / 44 / dz / dz 




From (19) follows 


t 


2 2 
•9*— _ n 


Z 1 " + w3. 
Z 


. R 2 - dw _ _ It 

t V * ; * TT 


2 - 5 2 


dz 


w'» &£ = - si--.. * 2 ; 5 k'. 1 2 1 

d£ z 2 . a2 dz ( z 2 . *2)2 


and elementary integration gives 

M §£ dz = (e 2 - a2) 2 / dz 

d z d£ j/ z 2 (z 2 - a2 ) 

/C^dz-(R 2 - a2) j £ ~ a2 d Z * (E 2 - a 2 ) In z + 


a 2 V 2a z + ay 


ai 

2z £ 


A *£ dz » - (R 2 - ft2) A 2 ,,* a 2 

J d2 / «3 


r 


dz * - (B 2 - • ; 2 ) ( In z - £ — A r 

2z V ! 


(22a) 


dz 


'* * 4. - a (E 2 - a2) £ / - 

d * / ( e 2 - *2) 2 

o 


_ (R 2 - a 2 ) 2 ( Z_ 


a 2 \z 2 “ " 2 


+ J- In - z -~— - 

2a z + a 


A ^ = 2 (p- - -. a ) / - ? - £ - t .. a2 . 
J 4, y (,a . a 2)' 


dz 


. 2 (R 2 - a 2 ) z 


z 2 - a 2 
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with (22) and (22a) , (21) can be integrated, giving as particular 
integral of the nonhomogeneous equation (21) 


4 W 2 * = A 1 + B v' 2 v + C w' ( r2 -&L .( \ + In f 


JL* a 


z 2a z + a 


+ (a + 2K) < a + *r *) 

+ (c * 4) v + V’ X b L J ^ (| + fe 


+ a 


;2 


dz d£ 2 


dz d£ 


v' 2 jK'AS JL* + ( ff + 4) w'M' is (^ 2 - a 2 ) ^ln z + 


a 

2z 2 > 


+ »'®3& I * § g' ag -~-a^ a -(| * k ln tffV * ,g w ( ° + * 

^ \ / 


- (a + 4) 


- 0 v 


,2 J^. 2 v i...^ll - (a + 4) y».CB 2 ..: A. / I 


-fi 2 
z - a 




2 2 
e: 


A_ l n £^Jl\ + $ *' 

2a z + a J 3 6 2 


dz d£' 


2z 2 J 


- 4 •, is v - (ff + 4) V« --.£)! 

d 4 .2 



( 23 ) 


The two terms multiplied by In z counteract each other, as should 
be, as otherwise an impossible circulation could appear in this 
approximation. For large values of fzf ** E, W z * behaves as 

VlO 
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It is advisable to rearrange this expression: 


VIq* * 


w' 3 j + «2 V 2 v + «3 v' + I; in ff&j 

' + h ^ tfi) + “5 v ' 2 w " f * a 6 v ’ w " (f) 


+ % ^ 


+ Srjr f ?’ 2 V + 


((_ r 2 


4ft 


*6 v ’ ST 1 1 + =9 w ' 2 ^“2 + *io *3 

z - a 


+ \i v w + *12 d£ w + ®13 




(24) 


vhere the coefficients to ta^ follow in simple manner from 

the coefficients of ( 23 ). 


In order that Wg comply with the boundary condition Wjg = 0 

at the image-plane circle, suitable analytical functions of z and 
2 must be added to W 2 *, as they, are homogeneous solutions of 
(21). For example, the term w‘~£ must be replaced by 
■ -il- b£ a 2 z t . a similar procedure must be used for the 

V V -* ’k 2 )’ 

other terms. Certain terms are accompanied by new singularities 
in the outside domain, which must be compensated for by adding 
further suitable homogeneous solutions complying with the 
boundary conditions- Thus, the coefficient of ao must, in 
order to satisfy the boundary conditions (at the image -plane circle 
and at infinity) be substituted by 





In 


z - a 
2 + a 



) 


and in addition, for removing the singularities at 
the following function has to be added 







-,2 _ 2 

R + ar 
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g2 _ ^2 

Since o o * a term ~ In e occurs) the same is the case 
I r + ar 

for a^. If one transforms all terms accordingly and makes a few 

simplications, a solution of (21) is obtained which satisfies all 
boundary conditions and has no singularities in the outside domain- 
This solution is Vq} its imaginary part represents the third 
correction tyj of the stream function: 
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as check it is easily seen that the residue of Wg of z, 

■n2 

z x. t disappears, hence the singularity at these points in 
(25) is only apparent; furthermore, in each term with. z z ■ r 1 2 
the factor E 2 - r 2 can he removed, henco Wg vanishes for 

j z | = J z | * EJ finally that the expression vanishes to a sufficient 
degroe at infinity. Computing the normal derivative of \j/g on 

tho ellipse from (25) and developing especially for z = i E, 
z =» - i B (the end point of the small axis) in powers of e , it 
is apparent that: 

(1) The development of all terms starts with a term of the 

third order in c in particular, at and through the 

logarithm, a term ~e3 (arc cos e - l) or in third order 

~C~ ^ - l^ will- occur that cannot he removed again h y the terms 
eq.ueQ.ly proportional to arising from \)/^ and \|/g. 

(2) In eg and ^ an additive term proportional to ^ In £ 
occurs likewise for which the same is true as in (l) . Also this 
term must enter in the final formula for the velocity in the 

end point of the small axis. The order of magnitude of this term 
is still sufficiently large for a small hut not too small e 
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compared with terms of the order thus it cannot be neglected 
for this ordor of approximation. 

It has thus been proved with this special exetople that within 
the known terms , the Rayleigh approximation and the Prundcl-Bnaemann 
approximation give the same result. 


Translated by J . Yanier 
National Advisory Committee 
for Aeronautics 
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PA«T II. EFFECT OF COMPRESSIBILITY 
IK AXIALLY SYMMETRICAL FLOW 
AROUND AN ELLIPSOID 


The differential equation for an axially symmetrical irrotational 
flow of a compressible gas is solved on the basis of the solution 
of Prandtl's linearized equation as first approximation for the 
case of axially symmetrical flow around an elongated ellipsoid) 
terms of higher than the second order in the interference velocities 
are neglected- The maximum velocities of ellipsoids of various 
thickness are calculated in dependence of the Mach number and 
the results are compared with those of Gothert (very slender 
ellipsoids) reference 1) and Lamia (sphere) reference 2). 


x, r, cp 
e* - d*/t 
U* 

0 • x,T * r/n 

e = d/t 
U * U* 

P, a 

P o > 

P i a 

CO CO 

v * Po/Pco 

M* - U*/a* 

u - y/ i - Hf 


S1MB0LS 

cylindrical coordinates in compressible 
flow 

thickness ratio of body in the x-, r- 
plane of flow 

free -stream velocity in direction of 
x-axis 

coordinates in the incompressible 
comparative flow 

thickness ratio of body in the c-, t - plane 

free-strcam velocity in the c-, t - plane 

local density and sonic velocity of 
compressible flow 

density and sonic velocity of air at rest 

density and sonic velocity at U* 

density ratio 

Mach number 


distortion factor 
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u, v 


component of total velocity of compres- 
sible flow in x and r direction 


= u - U*, v^ » v 
0>\|/ 

4>1, vj/ x 

it »i 
k 

T, q » l//l"-"c 2 


disturbance velocity 

flow potential and Stokes 1 stream 
function 

potential and stream function of 
disturbance flow uj_ # v^ 

potential and stream function of the 
free stream 

elliptic coordinates in the a-, r -plane 

scale factor of the elliptic coordinates 

elliptic coordinate of the ellipsoid with 
thickness ratio e 


C 




^0 + 3 - 
^0 " 1 


factor of the disturbance function 


K * c p/ c v ratio of specific heats 

I. RESULTS AND NOTES FOE THE APPLICATION OF THE PRANDTL RULE 
TO ELLIPSOIDS OF REVOLUTION 


In order to describe the flow around an ellipsoid of revolution 
cylindrical coordinates x, r, tp are introduced. The x-axis is 
placed in the free -stream direction, r is the distance from the 
x-axis, and cp the angle of rotation. The flow yields the same 
pattern for each <p * const; thus the differential equation of the 
flow is independent of <p ~ 
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0 is the potential defined "by irrotationality, u and v are 
the velocity components in the x and r directions, and \J/ is 
the Stokes stream function defined by the equation of continuity. 
If one substitutes the disturbance velocities 


= u - U*, v^_ = v 


for \j/ instead of the total velocities in the equation, where 
U* is the free stream in compressible flow, and expresses the 
local sonic velocity a by the velocities and the sonic velocity 
of the flow, one obtains similar to Busemann (reference 3) 
the disturbance stream function = \}/ - ^ for the two-dimensional 

case 


^lxx + ^Irr “ ^ (K + 1) ^ 


Vlr 

Ixx r 


05 


+ u 2 (k - 1) ~ 

" 2 ^lxr r T .2 
- (k - l) iu 2 


lrr r 

r 
C 

r 

*lr 2 


r 2 j 


( 2 ) 


M_ = — is the free -stream. Mach number, 
00 6 • 
00 

at rest, (p * P Q ) , and for free stream 
the stream function of the free stream 


o = 
(p 



S 


the density ratio 
P )» If ^ is 

03 CO 


U* a - u £ 
r 


(3) 


omitting the second power terms in the derivatives of the stre am 
function on the right-hand side, there remains the linearized 
differential equation of the axially symmetrical flew whose solution 
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corresponds to the Prandtl app roximat i on . By distorting the x, r 
plane with the factor p » l /\/ 1 - 

x = a r = pr (V) 


the linearized equation gives the differential equation of the 
incompressible ellipsoid flow, the solution of which is immediately 
indicated. Substituting the solution of the homogeneous 

equation in (2) on the right-hand side, the solution of the now 
inhomogeneous equation will give a solution ^( 2 ) and 

^1 + ^ i ^ (5) 


represents a higher approximation for the solution of the compressible 
fLow differential equation. The solution of the linearized 
differential equation, that is, the first approximation for the 
compressible flow, is in elliptical coordinates (Lamb, reference 4) 


c * k£n 


T =° k\ 


- 1 1 / 


n 2 - i 


.( 6 ) 


for an ellipsoid of thickness ratio e* and free -stream velocity 

U* in the distorted a-, T -plane € = — 

\J/ 2 . . 

U » - u*^ 1 ) 


C k 2 (1 - t 2 ) (r] 2 - 1) In . --Jl_ 


with 


-U/2 




2 ’lo " 1 


r 

(7) 


For i| b t) s —-==== the coordinates a, T represent tha ellipsoid 
\/l - € 2 

of thickness ratio £ and the first approximation of the velocity 
u-j_ follows from the solution \J/^ ip 
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r t Tj 


(8) 


The higher terms in the derivatives of the stream function are 
neglected. If one includes the square terms for the calculation 
of the second approximation, the first approximation 

\J/ gives an additional term A u^ 1 ^ in the velocity uj. 

of higher order 

u^(^) = + ^ (9) 


Substituting (7) in (8) gives the Prandtl approximation for the 
maximum speed (£ = 0) at the contour of the ellipsoid (tj « ti q ) 


U 1P 


(li 


max 


U* u 


2 2 


1 . + 1 1 


2-2, 2 T) 0 1 


( 10 ) 


with 


1 ~ ± 

l/l - (€*/n) 2 ~ 2 


Neglecting the terms of higher order in € for thin ellipsoids 

(e 2 ^ i) gives (asymptotic values for rj 0 -#l, first approximation^ 


L 


Uip 


( 1 ). 


U* 


max 




3 e* 2 fin 


( 11 ) 


With the same omission in incompressible flow 


C 1 H ‘t*(p *s-j) 

s u / max 6 


^•This result of the first approximation can be entered at once 
according to the Prandtl rule, when Gothert 1 s form (l) is applied: 
determine the interference velocity in the incompressible flow at a 
contour distorted with l/p and multiplied speed by (i 2 . 
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The ratio of the compressible to the incompressible maximum velocity 
for an ellipsoid of small thickness ratio € thus is 



In - 1 , In ji * 

£ * 1 + 

In | - 1 in | - i 


( 12 ) 


Consequently) the ratio of the velocities for thin ellipsoids is 
not unity in first approximation. Rather an additional term 
proportional to In |i is obtained, which at higher Mach numbers 
can be quite considerable even for small thicknesses- The value 
of unity is attained in the limiting process to vanishing thick- 
ness at fixed Mach number. Considering that the first approxi- 
mation of the potential by (7) is, like the velocity in 

(ll), proportional to e 2 , it is apparent from (2) that the 
additive stream function ^(2) can contain only terms propor- 
tional to e 1 *. Thus (ll) certainly contains terms all proportional 
2 

to € . Eor the second approximation of the maximum velocity at 
the ellipsoid with 



*o - 1 


- Qo 


as abbreviation, one obtains the following result 

U W ’ + u l (£) mM 

with 


u- 


( 1 ) 


max 


u- 


IP 


( 1 ) 


max 


A u, ^ 

1 max 


(13) 


Gothert obtained a similar rosult In p/ln A, as is 

apparent from the above when (in 2) - 1 is neglected with respect 
to In €• But such an emission is admissible only for oxtremely 
small £, because In e » - 2.996 is itself no longer great 
with respect to (in 2) - 1 = 0-307, so that errors of 10 percent 
result. 
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where 


|u (D 

IP 


max 


i 2 U* 


1o g - 1 

tn 0 2 - i) dc 


denotes Prandtl'e first approximation, while 


A ujfc) 


max 


H 1 * M 2 U* 

CO 


1 + (l + KM C 


2 ', 


(--*) 


is)- 

" (n 0 2 " ^ ^o 




and 


»l <a) Ba* ’ S 


1 M 2 U* 


2 - 


i 


- (% 2 " 1) Q 

1 + H 2 (1 + Mt®) 


n 

ti^ 2 ) 


~ 


n 

n=n 0 

- 


Tl^oi (HO 


represent the additional portions of the second approximations- 

H'XDr^ 2 ! 


The expressions 


!•»{; (2. 

and 

b m J 

JC!lo 


are solely 




dependent upon rj Q = l/i/l - {&/\i) £ t that is, on the ratio e*/ii. 
They are shown in figure 1 for email t) ^ 1 and therefore immedi- 
ately available for every practical example • t) 0 is the larger 

the thicker the ellipsoid and the smaller the Mach number. Even 
at e* ** 0*5 and Mq, = 0-3 tj is not larger than 1.138, hence 
still a number close to rj 0 = 1. 

Figure 2 shows the maximum velocity increments of the second 


approximation (u^/tr*) 


max 


U-- 

V u* ! max 


on ellipsoids of various 


thickness ratios c* plotted against the Mach number. The range 
of validity of the curves is bounded by the curve of the critical 
Mach number, on which the local velocity exactly reaches sonic 
velocity. The heavy curves of the second approximation 

are compared with the first approximation [u i: p^/u*| ^ aY , which 


max 
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follows from the first approximation of the stream function 

when neglecting the square terms in the derivatives of • 

This first approximation is represented by thin curves. To 
illustrate the magnitude of the two additional terms of higher order 
the total contribution by the linearized differential equation 

j U tl ^ ^ Aj*j max is included in dashed lines. It is seen that the 
additive term Au^^ contributes entirely too much and is 
reduced again to the greater part by the term of the second approxi- 


mation 


u-^2) , 


In order to Judge the quality of the second approximation, 
the maximum incremental velocities for the extreme case of the 
sphere d*/t » e# = 1 are plotted in figure 3 and compared with 
Lamia's values for the sphere at several Mach numbers. Since 
Lamia referred the Mach number to the critical velocity -of sound 
a*, the Mach numbers 0 = U*/r* were converted to the sonic 
velocity a<» at free -stream velocity U* by the formula 



y* 2 j _, 2gf 

e _ 2 K fc + 1) - (K - l)p 2 


It is seen that Lamia's values computed in fourth approximation 
by the Rayleigh method do not differ very much from the second 
approximation computed here; hence it may be concluded that this 
second approximation is surely sufficient for all practically 
encountered ellipsoids with thickness ratio up to about £* « 0-5* 

Figure 4 represents the conditions for the incremental 
velocities in compressible and incompressible flow around bodies 
of the same contour plotted against the Mach number for various 
thickness ratios of the ellipsoid in first and second approximation. 
Lamia's values are included for comparison in dash-dotted lines- 
Figure 6 contains the corresponding conditions for the total 
velocities u = u^_ + U* in second approximation* Figures 4 and 6 
represent the final result of the higher approximation which goes 
beyond the Prandtl rule. The compressible disturbance and total 
velocities of the second approximations are directly obtainable, 
if the incompressible values are taken from figure 8. 

In the calculation of the second approximation of the stream 

function and its derivatives the terms of the solution \|/, (2) 

IT) 
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vere represented "by infinite series which, for tj 0 — ? 1, tend 

toward logarithmic l imi ting curves- If one replaces these series 
for thin ellipsoids by their limiting curves and considers the 
terms up to the order one obtains for the incremental 

velocity the asymptotic formula (second approximation) 



= e *® (in f£ - l) 


+ [2 In 2 ft + - 3) In gi - | n 2 + 3 

♦ St (f -.1) ( i t (15) 

The asymptotic values for the partial solutions ^ and ^iin (2) 

employed here are indicated in figure 1. For the incompressible 
case for n =* 1 the asymtotic formula is true 



Thus the ratio of interference velocities is 




NACA TM No. 1233 


k2 


For the total velocity u = U]_ + U* the ratio is 


u 


SagE^—’ l a 1 + t* 2 ta w . ^ (n 2 - 1) 

^jncompr./ ^^ 2 \T 


- 2 ln^ln~ + 3 (in ^ - l) + |- <1 + K it 8 ) 

+ In n ^2 in ^ + %i 2 - 3^ (l8) 


The pimaM crosses indicate the asymptotic values of the second 
approximation for €* = 0*2 in figure 2 ; one con see that they 
are still a little below the first approximation at small Mach 
numbers, between the first and second approximation at medium 
Me, and then above the second approximation. Then they finally 

approach the second approximation asymptotically at ^-> 1 , 

since for p — >co in n * — ==*= 5 = 3 ======- and thus rj 1 . 

0 \/l - <e*/n > 2 

Figure 5 contains, aside from the exact second approximation, 
the values of the second approximation computed by the asymptotic 
formula (l) for €* » 0.2. Thus the two second approximations 
do not differ very much for small thicknesses. If one neglecto 
the terms of higher order in ( 17 ) one obtains the asymptotic 
values (ll) of the Prandtl approximation corresponding to the 
exact values of the first approximation. One can see from 
figure 5 that the exact values are too small, the asymptotic 
values too high, hence their average is close to the second 
approximation. The dotted curve in figure 6 represents the 
asymptotic values of the first approximation for e* = 0*2 
corresponding to ( 17 ) with neglect of the terms with c* 1 *} 
the dashed curve indicates the exact first solution. 
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II. OUTLINE OF THE THEORETICAL DERIVATION OF THE RESULTS 
(a) The Differential Equation for Potential and Stream Function 


The continuity equation of a compressible medium for axially 
symmetrical flow with local density p and the local velocity 
vector v is 

div (pw) = 0 (19) 

Expressed in cylindrical coordinates x, r, cp with v = (u, v, w) 
it he comes 



In the axially symmetric case the flew in every plane cp = constant 
is the same, hence only the flow in one x, r plane needs to he 
analyzed. In this plane (w s 0) 


|- (pru) + ^ (prv) = 0 


( 21 ) 


Thus there exists a Stokes stream function \| / with the property 


i P ° , 
tt ” r 


T-lS&y 

r p Y x 


If the flow is irrotational, there exists in addition a flow 
potential £> , hence 


U s $ s li 3^2 \l/ 

* r p * 


v =0. 



( 22 ) 


Thus from the disappearance of the rotation, rot = 0, there 
follows for the stream function 


ox \r p x 


+ d / 1 £o 

or \r p 




m 0 


(23) 



NAGA TM No . 1233 


Moreover, "by Bemouilli's equation 


+ ; q ^ t » Constant 
P 2 


with 


j 2 = u 2 + v 2 


.J 2 ) Js(+r 2+ +x ) 


and by definition of the sonic velocity r 2 *= d p/d p the 
relation follows 


^2 a d In p « - dq 2 
P 2a2 


Substituting (24) in the basic equation (23) gives 
s(^) + S?(!r) + r (+x + +r 


i-ia.p . - _i_ a [/" £a^ * *r 

&C 2a 2 « VP / y2 


1. /£a\ 2 tfxx £Mrx f d ln P <£ 

a 2 W a 2 


dlnp = 

dr 2 a 2 «* V P y 2 


1 /Pof ^ar+^rV + A . 5 In n a 2 

- 2 VP A r 2 ‘ ,3 / Sr a 2 
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or tetter 







The basic oquation thus is 






The differential equation for the potential is obtained, in the Bsme 
way. From the equation of continuity (21) and. the Bernoulli 
equation (24) there follows 


h (r ^ + b ( r 5 r) + r — fe- 


f r v . 




with 
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Hence the differential equation for the potential is 

$ 

®xx + ^rr + 2 '“’r ®rx + ^xx) 

a. 

- £ 2 4 > +<J } 5 A = 0 

fe 2\ r rr r x *rx / 

or, shorter 

(28) 

(h) Introduction of the Disturbance Velocities and Transformation 
to Elliptic Coordinates in the Plane Distorted 
According to Prandtl 

From the adiabatic equation p/p Q = (p/Po) t K “ 1*^5 (air) 
and the definition of the sonic velocity a -2 »• dp/dp there 

p 

follows when the pressure is expressed by q according to the 

Bernoulli equation, a 2 » a Q 2 - q 2 , Sq = sonic velocity of 

the gas at rest. 

By means of this relation the density can be expressed by 
the disturbance velocities 


Uj_ »■ u - U*j Vjl 


a V 


If a^, is the velocity of sound at free- stream velocity, q = U* 
and Pa, the corresponding density, there follows with ~ = u 

P« 


£_ 

Po 



K“1 




1 


L. 



p /, 

,2 \K-1 

ss 


) 


Po V< 

* 2 / 

CO 

(k 

- 1) 

U* U1 1 K - 1 

„ 2 2 2 



«co a 

00 


1 

K-l 
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4? 


or 




, U* u x x 

1 + 7T + 

a 2 2a 2 


OD 


03 


(l + K mJ ) Uj 2 + V-, 2 


with 


Jg2 


(29) 




2 


Splitting off the contribution of the free-stream velocity in the 
stream function there remains 


❖ = ^ co + ^l with 'l' oa “ - ^ U* r 2 


because 


U* » 



" u 


\j/» r 

r 


Now the disturbance velocities can be expressed by the disturbance 
stream function 

From 


* "l ■ ' f 2 r * ’t'lr) =r “1 " D * (i 5® ' X ) 

. 2s fc . M 2 _S2L Ri -K K 2 )u,a + T. 2 l -£a 

pr ^ oo o_ 2 I co 1 1 | p 


2 1 

follows with u = 


2a 2 L 
00 


£o ^lr 
P r 


1 - Moo 


% — 2 (l * + [(1 - « O »i 2 ♦ * 


( 30 ) 
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Hence In first approximation 


u^ 2 -i\t 2 


tiE 

r 


v, » 


„£i* 

r 


(30) 


If one substitutes this approximation (30) in the terms of higher 
order one obtains as the second approximation 


p^lr/, uu 2 U*^lr 
U 1 * r 1 1 ’ ^2 

V a « r ' 


2 

2&co 


(3D 


^ 00 


Substituting these disturbance velocities in the equation of the 
stream function ( 27 ), and neglecting all terms of higher than the 
second order one obtains, with 


^ r = -^U«r+^ lr ♦ IT -“2* + ^ lrr 


and 


a 2 st a 2 - (it - 1) U* 


U 1 


*lxx 


a 2 - (k - 1) U* u x - U* 2 - 2U* 


U 1 


+ * 


Irr 


hp 2 - ( K - 1) U* ui 


^♦ixr 


♦ 


lr 

r 


- (K - 1) U* u x 


U* o 

+ r T i - 0 
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1 “ M 2 + (K + l) tf|i 2 1 + (K 


U 0O «5 


ir !*(„.!) ffiL 0ll 8*l£ . 0 

r L r _ 


Thus under consideration of the terms up to the second order the 
differential equation of the disturbance stream function is 


(1 " M co 2) + ^lrr “ ~ u 2 ( K + 1) t 


$lE 
lxx r 


+ Jl 2 (K - 1) 


- a}/ - (k - l) u 


lrr r 


2^lr 2 


lxr r 


Distorting the coordinates according to Prandtl x = a r = ht, one 
obtains, because of 


*lr - j *t>* +irr ~ ^ Vr* *lx = V 


♦mr + ^tt " r 


* r 

I- 2 (« - U* ao t T + (" - 

t i 


*• * L 


^CT^C + — - ( R - 1) 


>1/(1) 


If V' ' is a solution of the homogeneous differential equation, 
that is. 


(1) Cl) Mv ; 

*„ 1 *+ T T - - 0 
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the Insertion of this solution as first approximation in the 
right-hand, side of (33) gives, with 


* 




for the correction to Prandtl's solution the inhomogeneous 
differential equation 


v CO T T 


(2) tr (2) 


. affiL 
2 


8 T 

CO 


(1) CD 


ft" 




U> 


- 2 * 


n 2 + i) 

( 1 ) 


at 


(K - 1)| 
(34) 


Since the solution (?) of the homogeneous equation iB given in 
elliptic coordinates, (34) must be transformed to the new coordinates. 
The simplest way iB achieved by not using the transformation (6) 
but rather by mapping the a, T plane at first on an auxiliary 
plane 0 by means of the analytic function 


g + IT a h cosh (£ + 10) 


(35) 


The transition to elliptic coordinates is then given by 


tj * cosh £ £ = cos © 


(36) 


The existence of the Cauchy-Eiemann differential equations for the 
mapping (35) facilitates the transformation considerably and 
(seo also appendix, section a) the traeeformaticn by means of (36y 
involves merely rewriting of- the result with the new notation. 

The differential equation (34) attains then the form 
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<i 2 - x) * m < 2 ) + (i - 5 2 ) * u <2) 


2ES 2 — — (1 - tsj 1 
•2 (n 2 - 5 2 ) 2 1 4 


V. 2 (IS + 1) - (« - 1) 


+ 3 ( n 2 . D - ^(n 2 . -J) 


,2.^2 


2 . 


(*l 2 - t 2 ) n In ^-7^ * 2r, 2 


+ 2g‘ 


-^ + 2 (^4t) 2 - £ 2 - 1 - z p\ 
1 \ ^-S 2 / 


(c) Solution of the Differential Equation 

In order to solve (37), the equation is "broken into tvo 
differential equations. 

(„ 2 - 1) ♦ (i - S 2 > ^ - *£7^5 [t* - 5 EV 

+ 3?a + ” 2 ] - 11 ^ (2> a 


+ 2(n 2 - x) - 5SV + 3n s + E 2 ) 


and determines one particular colution for each of the two. 
is accomplished "by solving the equation for individual terms 
the right-hand side until finally all terms on the right are 
exhausted. In this way one obtains as particular integrals, 
the simplification (in x) 2 * ln 2 x 


£ 2 


(37) 


(38) 

This 

on 

with 
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This solution substituted in ( 38 ) gives exactly the terms on the 
right-hand side- The indicated solution vanishes as l/i-p aside 
from, a constant at infinity, hence at infinity parallel flow is 

attained- To satisfy the boundary condition \p 2 ^ * 0 at the 
contour tj = t) q , solutions of the homogeneous differential 


equations must be determined which added to the nonhomogeneous 
solution makes the latter equal to zero* that is, cancels the 
previously obtained terms (39) for t} = tj 0 . This is accomplished 

without difficulty for the terms of tho form (l - £ 2 ) cp(rj) and 
, now we set up the equation 



* (l - t 2 ) 



(n 2 -t 2 ) n 


If the differential operator at the left-hand side of ( 38 ) is to 
disappear, the condition 
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- 1) <p" " 2(n - l)(2n - l)<p 



r 


1 

+ 4n(t] 2 - t 2 ) ' 

- n(n 2 - 1 ) q>’ + 

(2n - 1) rj 2 - n 

r 


m 


j 


will exist- For n ® 0 the differential equation is 


a 0 


(n 2 - 1 ) Cp" “-2<P a 0 


that is. 


<P(ti) = c x | (tj 2 - 1) In - n 


+ Co Cn 2 - 1) 


If n « 1, "both brackets must disappear, since the solution 
<+>(t}) is to be independent of £ - The oquations <P" =* 0 and 
rj cp* - cp =0 are both satisfied by <P ** C^- Hence one obtains 

as homogeneous basic solutions: 


- a - £ s ) (-1 2 - 1), % - a - : e > 


mH4-i 


(40) 


and 


if, - (1 - £ 2 ) 


i s - i 2 


But for n = 2 two differential equations are obtainod which have 
no common solution. Thus for tho terms of the type 


1 - i 2 

(n 2 - £ 2 ) 2 


q?(n) 


and 




i 2 
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no closed solution complying with the homogeneous equation can 
he given, because the second expression, also which can he 
reduced to the first by differentation, is of the same type* 

These expressions must therefore be developed in terms of 
functions which are solutions of the homogeneous equations- 
The general solution of the homogeneous equation which disappears 
for 11 — is of tho form 


* 


I 


n 


n(n + 1) 


(i . ( S) („a . i) jJSeM 


where P n (£ ) and Q^n) represent the spherical functions of 

the first a n d second type. Thus the above expressions must be 
developed in series of this type for t) ** rj . A homogeneous 

i n 0 2 * 1 

solution which for tj * t] 0 assumes the value - p In — r r 

v " L 

is obtained (appendix, section b) by setting^W- «= P «'(£)* 

- w 1 

. „V~(D 4n + 3 O&u-lM 

* 4 "vSS (2n * u & * £) (lo a - tio 2n+1 M 


u - 5 2 ) 4»i <£) (I 2 - 1) 4*1 (i) 


Ctt) 


Since this series converges absolutely and uniformly for q £ t) q 
and (£(«£ 1 (see appendix, section b) for n =1 1 0 

> M Tl o) s " 2 ln “ 2 I £2 


Differentiating twice with respect to n gives a new series with 
the aid of which the second expression can be also represented 
(appendix > section b). 


X 


( 2 ) 



02n+lfao> 

^ ^ 2 n+i Cn 0 ) 


U ' t 2 )^n + l UK* 2 - 1) <4n + l M (4S) 



The series also converges uniformly and absolutely in the same range and two homogeneous solutions 
give 


te(, 0 ) , ; 3--- £ !£V <2) 
* 5 ' *<>' 2 2 .2 2 ^- 

n * c 


and hence for tj c ^ 


» n 0 3 —~ 0 " t 
* (n 0 2 -C 2 ) 2 


Having thus the terms of the inhomogeneous solution supplemented by homogeneous terms with 
corresponding constants so that the expressions for jj * rj 0 disappear, these expressions are 

then appropriately expanded in powers of i) - t) q . The first tern then gives the derivatives 
of the two parts of the stream function with respect to rj at tj = « , needed for computing 
the maximum velocities at the contour alone: 0 


a*j{g ) 




- (1 - &) 


lw 0 


= (i - £ 2 ) 


2 ln 2 1)0 + 1 . i 
_ % - 1 1 

|(i, 0 2-l) in^li-Tj/% 2 ^ 

- l o 

1 1 o ^0 "** 0 1 + 1 

4 1x1 Up - 1 ~ gT »o 111 rj 0 - 1 + ^ - 


o \ + 1 

1 x ^Q 3ln n7^i f^o.0 


2n, 


3SZJ 


(n 0 2 - S 2 )' 


3 + 4n c 2 + tTT^ i - jl 


i (^o 2 ~ 1) In *° 


+ l . 

- 1 ^ 


*io 2 - £ 2 


gFl(Hp,S) %_ 


% a - 1 2 (n 0 2 - 1 ) fa 0 2 - ? 2 ) 


+ a_ V (1) 

3ti j ~ Tj=rjf 


C^3) 


vs 

VJ) 
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with 


>1 (lo,f ) ■ 5 1o 3 - 3i 0 - 5 ^ (lo 2 - 1) (lo 2 - £ 2 ) - >H1 0 3 


V ~ 1 
1 2 - £ 2 


| (% 2 ■ £ 2 ) 2 


4. \~ (2) 

^1 ^^0 


(n D 2 - 1) 


and 


3 y (1) .£ (ta + 


ST 


^n+lflo) 

3) tt - £ 2 ) p ^ +1 (£)^ n+ 1 (n) 


(% 2 - D ^ n+1 (n 0 ) 


a r-( 2 ) y— 


^Z_ 


2 _C^n + 3)(2n + l)(2n + 2) 

, o 

Ci<, 2 - 1 ) <4»i (n 0 ) 1 ' £ ,P2 "- 1 U) <W ’ ,) 


The derivatives of the aeries follow from (4l) and (42) when 
observing the relation §- ^2 . l} V(n )j i n(n ^ 

(reference 4) . From the derivation of the Btream fraction ^( 2 ) 


W£) 

’on 


- »S? 2C 2 * 2 

a 

00 


i l 


u 2 (k + 1} - (k - 1 ) 


5^(2) 

arj 


a%( 2 ) 

-*n 


*^ UV 4 


(no 2 : l ) 2 


1 / n 2 

o n^o 


1) In 


n 0 + 1 


no ‘ 1 


there results the velocity 
following relation for the 


according to ( 31 ) 
ellipsoid center 


w 2 A®i 



by means of the 

(£ = 0) 

3\l/ 
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since from (6) follows 


0 - ktt T n + n T £) m £ t 


1 = k 


- 15 7i 


."by differentiation, hence with T 


(£“0) 


k/n 2 - 1 


— — ^ -==== (>jy f g, + Tyr)—*^ ^ 

T t«o kv 4^Ti V* k 2 n 

The correction u-j_^ 2 ^ for the velocity in the center of the 
ellipsoid therefore is in first approximation 


^ (2) . _ „ 4 


,• u* 


U*r 


% 2 - 1 


n 0 - (% 2 ; 1) 


^OfU-Ot-iT 


Substituting in its stead the first approximation in 

(31) for the velocities, one obtains 
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The total interference velocity in the center of the ellipsoid is 
thus 



•m s (1 + « m/)J 


+ tca<!> / 

T) l C 


Substituting in it the derivative of the stream function 

according to (43), it is seen that in the practical evaluation of 
the formulas only the calculation of the two infinite series 
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offers any difficulties. The derivatives of the spherical functions 
p n u> can be given at once for £ = 0, vhen bearing in mind that 


(£ 2 - 1) r n d) = n(£P n - P n _!) (reference 4) 


one obtains 


p S»l < 0> - < 2n + x) p 2n<°) 


(2n + 1) 


1x3x5. . . (2n - l) 
2 x 4 X 6 • • • 2n 


(~l) n 



(47) 


The spherical functions for small n are best computed by the 

recursion formula 

** 

n (^ « (2n - 1) n 0 Or * 1 “ (n - 1) V2 


with 


Q 

^o 


1 

2 


In 


+ 1 

3 " 1 
•o 



.fto * 1 


!> 


(48) 


vhilo the derivative is obtained from the difference of two 

Q n : 


(ti c 2 - 1) Q n ' - n (tj q Q n - Qn-l) ( 4 9) 

d * 8 V*( 2 ) 

The convergence of the Beries and grows worse 

as the ellipsoid becomes more slender and the Mach number greater, 
that is, the closer t>- 0 approaches unity. Since these values are 

of particular interest, a larger number of series terms must be 
calculated. However, even in the calculation of Qj^ 2 ^ to six 
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valid digits "by the recursion formula the initial value to 

20 digits must he exactly known, while, to obtain for 

three digits exact, Q n to is required. The higher ^ 

with n > 10 must therefore he determined by a well converging 
series* As the conventional representation of tho Qj^Ho) as 

power series in terns of — converges very slowly & new 

*)o 

representation as hyper-geometrical series is used, which 
converges particularly well far higher n* (See appendix, section c) 


The infinite series 


a y~U) 


and computed this 


way for a number of rj 0 values were approximated by curves, which 
pass well through the calculated values for the individual ■ T) 0 


and for ni 0 -•* 1 and T) q have the limiting values of the 

considered series (ipperJLix, section d, "(79)* The accurately v ccmputed 
values are given in table I, the approximation curves in 
figure 7* 


TABLE I 

calculated values for tee infinite series 


“ 

3 VW 

W-~ 

s r(2) 

"Stj 

q o 

i_ a_y< 2 > 
q Q 

1 

1.25 

1*5 

1*732 

2 

3 

-0.495 237 
-.219 56l 
-.126 722 
-.076 117 
-.020 096 

0 . 08 l 508 

-.125 57 8 
-.120 631 
-.093 211 
-.033 401 

-1 

-.450 7 84 
-.272 842 
-.192 447 
-.138 570 
-.057 984 

+4 

.074 192 
-.156 052 
-.183 197 

-• 1 6 Q 688 
-.096 375 


A closer analysis of the behavior of the infinite series in the 
vicinity of t| 0 » 1 yields as asymptotic representation for 

small tj 0 (si) (appendix, section i, (76) and ( 78 ) and £ = 0* 



( 50 ) 
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These approximations for rj valid near unity are also shown in 
figure T and it is seen that even small departures from tj = 1 
are accompanied "by perceptible differences in the values of the 
functions- An exact calculation therefore requires the use of 
the interpolation curves of the exact values even for small t) • 
But if only a rough estimate is roquired can he set equal 0 

to unity in the additional term of the second approximation 
except in Q q , and the asymptotic representation of tho 

derived sums 


a VW 



a r {2) 

di}L_ 


~ ^ o 


(50’) 


may "be applied- In this case the expression for the maximu m 
velocity at the ellipsoid can he simplified. Inserting (50') 
in (43) givos at £ = 0 , for the derivatives of the contribu- 
tions of the stream function the asymptotically valid limiting 
curves (fig- 1 ) 

^ Z) /o + = -«o 2 + - 5 (51) 

/ ' 

For near unity the following developments are valid 



Inserting these expressions in (46) one notices that only the 
term of the lowest order needs to be considered in the, parenthesis 

because the entire bracket iB already proportional to ( £—y. 

/ 

Therefore tj q may be Bet equal to 1 everywhere in the bracket 


with exception of the terms 


ln^j 

£* 


one thus obtains up to 
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+ 



2(n 2 - 1) In |± - | fa 2 - 1) 


+ ^ - fc2 - » (i + K m 2) 
2 00 


hence 



.*2 


max 




+ iJ* 

2u 2 


2 In 2 % + (2JJ 2 - 3) In % 


e* 


j£x 2 + 3 + ^ (l + K m 2 

c P 09 


( 52 ) 
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III. SUMMARY 


(a) The effect of compressibility on the axially sy mme trical 
flow around an ellipsoid was determined by computing the maximum 

TJ* 

velocities as function of the Mach number M » (a = velocity 

05 ^in 05 

of sound at free-streom velocity U*) in second approximation for 
ollipsoids of various thicknesses. The values obtained arc compared 
with the values of tho first approximation which arc obtained by 
computing the incompressible velocity values for a slender ellipsoid 

of thickness ratio r * \/l - M^ 2 ~ (&* *» thickness of ellipsoid 

in compressible flow) by the Prandtl-Gothert rule (reference l) 
and the result is multiplied by tho square of the factor 

U » ~ -i! opposed to the initial distortion. 

VI - M 2 , . 

For small thicknesses ( “ K 0.1) tho socond approximation 

almost coincides with the Prandtl approximation, and the 
differences ore still small up to sZ ^ O.J (error < 5 percent 

"fc 

of disturbance, velocity at = 0*8 (fig. 2). 

(b) A comparison with tho maximum velocities for the sphere, 
calculated by Lamia with the fourth approximation of the Rayleigh 
method indicates that even in this extreme case tho departures 
even of the second approximation are not appreciable (error 

** 2.5 percent of the disturbance velocity at M ** 0-5). Thus 
the second approximation is still amply satisfactory for thick 

ellipsoids (i2< 0.5) (fig* 3) • 

(c) The ratios of the interference velocities 

u lcompress./ "Uncompress. 011(1 of the tQtQl velocities 
"compress./ u incampress. 0X9 represented in diagrams (figs, k 
and 6) for the point of maximum thickness. The velocities of 
the second approximation calculated here can be obtained directly 
by these diagrams, if the incompressible values oro token from 
figure 8. 

(d) Developing the first and second approximation for the 
maximum velocities in powers of the thickness ratio (the logarithmic 

quantity In (^-) being regarded as of the order of magnitude of 

unity; gives expressions of first and second approximations 



64 


NACA TM No. 1233 


asymptotically valid for small thickness, depending upon whether 

/ d *\2 

only terms proportional to or also those proportional 

to (^~) ore included. For small thicknesses = 0 . 2 ) 

the second approximation con "be replaced hy the clearer 
asymptotic expression without introducing an appreciable error. 
In this manner simpler formulas for the velocity ratios of 
compressible to incompressible flow are obtained. In first 
approximation one obtains for the location of maximum thickness 


u lcanmres3. - 1 + In u 

u lincompress. ^n 2 1 

d* 

t 


and 


“■aWCMfc. « 1 + (i*) 2 Inn 
Uncompress. ' 


IV. APPENDIX 
Auxiliary Calculations 


(a) Transformation to elliptic coordinates ." The urans- 
f ormat 1 on (6) 


a e ay t = J3S 


with 


3 = k sinh I 7 * cos 9 6s sin 9 

p 2 = k 2 7 s ♦ 8 ^ . 1 


a«k cosh | 
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gives bocause of the Cauchy-Rieaann differential equations 

sT e “ °e * “ T g- -«»s °n • T el " "°ee “ c 

= - T ee * pB 

Conversely with 

D « = p^ + ]j 2§2 B a 2 _ ^^2 

s 0 -°,*r *-< T — ? 

£ o„ -&«£$, +£ cs ® 0 - ^ (» * 2» 2 ) I 2 + *§ (b - ^ £ ) y 

, -S 

X)2 y 3> / TO " > TT 

%T =£„ £ £ t - t„e»T - g (» - 2*®) g6 ♦ g <D - Za 2) 

•$(***- 

For the stream function, then 

= %% +^ e e a = | (M'g- otsj/g) > * 6 e T = A (as^ + py*l/ 0 ) 
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and 


era 


a 2 5 2 - B 2 ? 2 
D 


^aa 85 %s S a + ^6 e a ) S a +* s £ ff(J + ^e© e a) e o + V« 

* + ^66 a2B2 ' 2 % 0 a37& +^g«P ( 7 2 - S 2 + 2- 

-V 8 (f + p 2 - 

^(JT + ^0 e T )&cr + ^gSaT + (^egS T +^e© e T) e a + 'Mar 

= ^ p 1375 ^^ " ^ee) + %e^ 2y2 " ^s 2 ) +^7o(a- 2 + p 2 - iis£i5) 


>(53) 


(- 


+ i]/ 0 ap ( 7 2 - 6 2 + 2 ~ -- 2 - ^ - 



The left side of the differential equation (3b) thus becomes since 

+ v - -i (* H * * 9e ) 

W 2) - (2) +t ee (a> - ■ft5 (S) - 


(54) 


vhile 




( 1 ) _ 


KP 2 5 2 


U In SL±. 

\2 a - 


Sfe )K 
P 2 






2 . 


% " 1 
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must be inserted on the right-hand side of ( 34 ), so that 


apt ■ KP 2 6 2 


(a£ In - 2cfcJ * K7 2 5 2 (p 2 li^±_| - 2ak^ 


% 


(i! 


U 


K5 £ 


(a 2 + p 2 ) In SL±j£ - 4ak 

ft " 1C 


(1) 


t ee '~' * K(y 2 - B 2 ) ^p 2 in - 2ak) 

« K 76 ^2a01n - 40k) 

Inserted in (53) one obtains 


^ (1> = 


(a 2 + P 2 )7 2 + (7 2 - B 2 )a 2 - 4a 2 7 2 
+ a2(r 2 - 8 2 ) + 2a. 2 . 7 2 (a 2 + P 2) + . 4 <£££ 


- 2kaS £ 


2 p 27 2 + (y2 . 5 2 )a 2 - 4p 2 7 2 + P 2 (7 2 - 5 2 ) 


+ 2p 2 ft 2 §?. .f ft 2 ? 2 - 7 2 (a 2 + p 2 ) + 4 
■ \ 2ka5 2 (l> + 30 2 - itS^S 2 ^ 


y2 o2-,2 
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^ aP7&ln js 2 (a 2 + p 2 } - p 2 (7 2 - S 2 ) + 2(p 2 7 2 - a 2 & 2 ) 

+ g 2 (a 2 + ^ . h£p£ + p 2 (7 2 - 6 2 ) + 2P 2 

- 2kp7S jj2a 2 5 2 - a 2 (7 2 - 5 2 ) + 2(p 2 ? 2 - a 2 & 2 ) + 5 2 (a 2 + p 2 ) 

- + a 2( r 2 _ & 2) + 2a 2 s£sl^-&£ i 

- - ^2 ^P? 8 ^ * 2a 2 + p 2 - — 

. s ( D 3 ^ . aofcVts^. 25142 s 

T D\ a-k / T DP 


The thus obtained expressions can he immediately transformed into 
elliptic coordinates- There is: 


SL a r\ 2 - £ 2 ^ since S »• tj; 7 «(, and P. 2 *» . ( 2 - 1; & 2 = 1 - ! ' 


Now, since 


% ’*<m $ *\ f 

+«e ' V 


The left side of the differential equation "becomes 


fc 2 ^ 1 . £2) ^ ' X %J + 
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and the expressions on the right-hand side are 


4t1 2 (t 1 2 - ll 


,2 . £2 


| ( 1 2 - t 2 > 1" ‘ 2, l 




■ 2 -i 2 

i ) 2 -1 


+ 2 


2 - S 2 + 2 n 2 + n 2 - 1 - - l) 

n 2 - £ 2 


Thus the differential equation reads: 


u 2 - i» ot (2) + a - : 2 »f f (2) - - 2® 2 


A 2 


U 


v u 2 - i 2 > 


(1 - £ 2 ) 


V 


11 2 (k + l) - (k - 1)1 R 2 - t 2 + 3(tj 2 “ 1) 


_ 4n 2 (n 2 - l) 

ti 2 - i 2 J 


( n 2 - i 2 )n In 2L±1 - 2if 

T) - ± 


, 2 2 J t 2 - £ 2 , 2 L 2 . f 2 _ 1 _ iaSCflL^Jl 
^ [n 2 - 1 r 5 n 2 - £ 2 . 



(5 


(h) DeveloTiments in spher ical functions and _c onvergenco 
studies for the series .- To obtain the development of 
2 

- in 2 — I— =L- at ti = in a series of the form 
TJ 2 - L 2 




4 -: 2 ) p n '(£K» 2 - 1 ) 
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multiply 


% 2 ~ 1 
n 2 - 1 2 


’YjrS 1 ' 


with P n '(£) (where ^ *» Aj-^tJqS “ one then obtains, 

owing to 

’ 0 a m 

(1 - £ 2 )? n ' Pja' = < 2 2^ n + + 1 ^^' n = m (reference 5) 
for the series coefficients by partial integration 


2n + 1 
2n(n + l) 


I tt\ ^O 2 “ 1 


P„*(S) in 


2 _ £2 


’ + ^2 4£ ’ j W n^fi) » - *»' 

I* 

* (,) = i/.i p " a> r?S * I £ p - ( - £) -/ft 


n » 2n* 


As for 


\n+l / 7 1 


)1 e * 

p n (£) 

-i * 2 - J 2 


Hence the required homogeneous solution which has the desired 
value for t) «= tj is 
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V) 


-tx) 


03 


n«=0 


iffil + a 


^2n+l 


(2n f + l)(2n' + 2) (^2 . (tj 0 ) 


(l - £ 2 ) P^n'+i (^)(n 2 ’ ^Qsn'+i^ > “ ln 


r\^i c 


2 -l 

, 2 -S 2 


(57) 


To obtain 'the further development of (1 - £ 2 ) r-p, the 

,,2 ^ ' & 

homogeneous solution t) r~ ~ a — is represented. "by a series that 

T} 2 - £ 2 

must converge for all t} > 1: 


n 2 - £ 2 


y- =;'y(i - £ 2 )p n *a)(n 2 * DV^ 


For t] =* tj 0 follows 


Tl 0 (l - £ 2 ) — 


„ 2 _ f2 

^o ^ 


«) B n Ki -£ 2 )p n i a) 


with 


V - V(lc 2 - l>V(t!o) 


and as before, the coefficients follow "by partial integration as 

■d * 2n + 1 f 1 p 1 - £ 2 

B n " 2n(n + 1) n 0 ( P n ^ 2 _ £2 


2n + 1 f :x 

2n(n + l) T, o # / 

’•J "1 

/* 

0 


1 2£( no 2 - 1) 


C% 2 - £ 2 ) 


2 P n «) ^ 


n = 2n' 

= )’ ‘ ^fSrir (ll o 2 - l)Qn'(n 0 ) n = 2n' + 1 
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By (56), integration with respect to z gives 

11 p It) is. a £ J° * « ven 

-l n (z 2 - i 2 ) 2 i ‘ Ci n' (z) n otW 


(58) 


Hence for all i] > 1 the series development 

1 .3. 


no -t g ) _ \ 


2 . £2 (2n' + lK2n' + 2) 


(1 - £ 2 )P' 2n . +1 (0(ri 2 - 


(59) 


is valid. Since 


d_ 

dtl 


(t| 2 - l)Q* 


2n 


■*i] - 2 <“’ 


’ + + U^n'+l 


it follows "by differentiation which is permissible for ^ > i] 
that 

„2 " [2 ^2 _£2)2 " ’ XZ ^ n ' + 3 ^ 1 “ ^^n’+l^^n'+l^ 


( 60 ) 


For tj =* t] q there follows multiplication with ti 0 from (60) and 
subtraction of (59), 

_ )~1 Ur' * 3„ ( f 2 )r , ( f } 

( n 2 _ ^ 2)2 " '-g-S (2n' + l)(2n' + 2) ^ & ;P 2n+l U; 


%Q2n* + l(no)(2n' + D(2n f 2) - (^2 - DQ^,^) 



On the other hand the following development Is to he valid 


i 0 3U - t g )'i 

K 2 - ? 2 ) 2 


■ ? 2)p W + id>(n 0 2 - - C 2 )r n '(?) 

0 0 


(a) 


By comparison of the coefficients one therefore obtains owing to the Legendre differ 
ential equation (reference 5) 


tl o Q 2n , +l^o^ 2n ' + 1 H2n‘ + 2) 


A** 

2n'+l 


"1 


+ %(% 2 ~ 2n’+l^o^ 

1 4n» + 3 ^o^ + 2n l +ltnp) + ^c/'Ho 2 " 8n'4-l^ 

2 (2n ' + W**' + 2) (n 0 2 - i)« , 2n . +1 (\) 


o ) tm 


This development can also he secured directly hy starting from (6l) multiplying hy P n '(£) 
and integrating over £. The integral on the left-hand side 



- & 
fa 2 - t 2 ) 2 


p n *(e> « 


gives after partial integration 
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and these integrals can be oxpressed according to ( 58 ) by the 
derivative of z) . By differentiation for n ** 2n' +1, 

( 58 ) gives 


Q "(z) 
n ' 




kz‘ 


( z 2 - £ 2 )3 


a£ 


( 63 ) 


hence 




2n(n + 1) 
2n + 1 . 


2 n 0 2 V “ ^o 2 - 1 ^n' + '0 o (n o 2 " l)Qn" 


thus the relation (62) is obtained. Adding the serieB (59) and 
dividing the series (6l) by -2 for arbitrary t] one fi nally 
obtains 



2 i)2 - t 2 


) jtB. t ,3 
A— 2 




(Hq^ " 2n!-l^o^ 


a - S 2 )p’ 2n+1 (D (11 2 - i)Q' 0 _,(n) — * 


n 0 3 (l - £2) 


2n+l 


C% 2 " £ 2)J 


(ft-) 


Whether the obtained series 



and 


y-(2) 
/ 


actu a lly attain 


the values given in (57) and (ft) must be proved by absolute 
^nd^ uniform convergence. Differentiation by terms gives the 

DX ^- |jl ( ^ neoded for the velocities in (43), the 


convergence of vhich must also be proved. It is seen that 
the (l) series differ from the (2) series only by the fact that 



„ ( 2 ) 

(2n + l)(2n + 2) < 


is valid for the nth term ctj^. Therefore only the convergence 
for the (2) series needs to bo proved. To this end, the spherical 
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functions are to he estimated, first* From the integral representa- 
tion for Qj^z) (reference 5) there follows: 




1 fl (1 - t 2 ) c 

2 n+ l J (z - t) n+1 


„_i_ f 1 . a-t) D z-t 

S”* 1 ,/-! (* - t) n (z - t) 2 


(1 + t) n dt 


< 



i - t \ n z + l 
* 2 - (z - t) 2 


2 11 dt 


Putting 


1 - t 
z - t 


u 


du = z - 1 
dt ‘ (z - t) 2 


and 


( l - t\ n z + 1 
V Z " (z - t) 2 


dt 


z_±_l v n 
z - 1 ' 


du 


gives 


Q n (z) < 


1 z + 1 1 

2 z - 1 n + 1 



P4.Y 1 1 

V 2 + 1 / (n + l)(z - l) 


(65) 


and in addition 


(z 2 - l)Q n '(z) 


ajU, f 1 Ll - t g ) a 

211+1 J - 1 U - t> n 


dt 


( 66 ) 
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For, from the differential equation of the spherical functions 
there follows by integration 



l)Qn (%) r) 0 «sl' ’ ^ + 


/ T| 

0 ^(z) dz 


. pfa ±. , i) f 1 a - t g ) n 

2,1+1 J 1 J-x (« - t)"i 


dt d2 



hence 


(z 2 - l)<V(z) 


*Ltl f 1 fl - t 2 )" 

giw-l J (z ” t) ri+ ^" 


(z - t) dt 


and finally, since 


(1 - £ 2 ) 



Observing that z » 
that 


> (n + l)(z - l}Q n (z) (g T ) 

? n<^ j < 4 at { < 1 

^<2n (68) 

T) in the denominator in the integral (66), hence 


n ^ p n-l * £ ?n> < « 


3 

•n-lj 


tl* 


n! 





NACA TM No. 1233 


77 


it is seen to tie sufficient to prove the convergence of i ^ 

V " (2) 

for tj c tj o . The same holds true for the series r- ) , a a 

from the integral representation of Q n there follows corres- 
pondingly 


^2n+l^o^ 


< 1 for > t) 

o 


Therefore, only the series 


)J 2) n=ti 0 - )_/ 4n + 3)Q2 n+1 (fto)(l - 


and 


a v-(2) 




} (2n + l)(2n +• 2)(4n + 3) 


^ 2 2n + l^ 0 ) 


(fto 2 " U)Q'2an(n 0 ) 


- ^Wl (t) -L** 


needs to he analyzed. From (65) and (68) there follows with 

tJ^Vi * * < 1 for ft 0 > 1 



< iajL. 2 a 2n+i 

n Q ' 1 * 


(2n + 2) 


2(2n + l)< 


- — -W (2n + 2)g Sn+1 

Mo 


Hence \ ^ 
|Cl < 1 and 


is thus uniformly 
ft £ ft Q > because 


and absolutely convergent for 
a converging comparative series 


with greater members can be given. Corresnonclinglv owing to 
(67) and (68) " 
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1AJ < XSEL± + 2 ) n („)(i.f2 )p i 
1^1 < (2n + 2)(n 0 - 1) ^n+l^oKl l >* 


< (2n + l)( 4 n + 3) 


-1 


2n+l 


(2n + 2)( no - l) 2 


2(2n + 1) 


Hence 


M < (t| ^ ~ l) 2 ' (2n + 3)(2n + 2) ^ 2n+1 


The comparative series however converges, with certainty, as second 
derivative of a power series for q. < 1; thus the absolute and 

jS /g^ 

uniform convergence of ^ ' itself is proved. 


(°) Calculation of the spherical functions Q^z) for the 
argum ents . gj <* 1 and n 10 by means of a hvoerpeometric 


series . - For z > 1 the integral representation 



x/z 2 " 


1 cosh d$ 


( 69 ) 


is applicable. Introduction of a new variable X by 

(z + i/z 2 - 1 cosh ^ (z - \/z 2 - 1 cosh x) * 1 

so that 


ax 


= Z + 



1 cosh & 


gives 



\/z 2 - 1 cosh X) n dX with X Q « | In j- 
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Through the substitution 


h » z - \Jz 2 - 1 cosh X 


Qn’( z ) * I 


iA 


L n 


dh 


J 0 \/l - 2hz + h 2 J ' 0 


h n (l - hx)~ 5 2 dh 


with 




2 1 

X * 2 + \/z - 1 “ « Z 


-/ 


z 2 - 1 


Putting h = ^(l - o>) the integral representation of a hyper- 
geometric series reads 

la-^(x- r ^_ ’ 


2 x It- . . .2n 


3X5- • -2n+ 1 


(X 2 - 1} 

/ 

„ . 1 1 1 
F : n + - 


t 2'2 


x/l 7 


where F indicates a hypergecmetric series with the variable 
X » Z 

F fa,J3,y 5 — ===Vy a — » 1 + — r — 

\ y/l “ * 2 / V ^ (1 - z 2 )^ 7 X 1* 1 - ^ 

. ftfe+iMlL+iJ- k_ + 

7(7 + 1)2-' 1 - z 2 “ 


g(a + l)(a + n)S(3 ± l) (Q + n) 1 

7(7 + 1) ... (7 + n) (n + l) (1 > z 2 )n+l 


where ^(t^) Is replaced in the series by z = r) Q + \f t ) q 2 - 1. 

It converges for z 2 - 1 > 0, t] q 2 > 9/8. This series converges 
very quickly for greater n and not too small tj q . The QjJho) 
with n > 10 are readily computable by this formula, for 
n < 10 the recursion is proferable. 
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The coefficients of the series 

(l X 8 • . .fgfu - I) * it] 2 


<v 


( n + i) ( n + i) ■ * * ( n + gJ i i ~ ) s2ii fd/ 


(72) 


are computed by recursion for n = 1 , 2 . . -19 and 
H = 1 . . . 7 • There is 


^i+i ® 


with 




(2u + l) g 

( 2 (i + 2 )( 2 fi + 2 + 2 n + l) 


The obtained numerical values are given in Table 2* along 
with the factors of the series for different n 


for n 


An * 2 


2 x 4 ■ , , 2 n , 

3 x 5 . . .. 2 n + 1 ’ ^-+1 


2(n±lJ 
2n + 3 


A-, 


1 . . . 19 (table 3)- 


(73) 


TABLE 3 


FACTORS FOR THE SPHERICAL FUNCTIONS ^(z) 


n 

An 

n 

An 

1 

1.333 333 333 

11 

0.517 019 48l 

2 

1.066 666 667 

12 

.496 338 702 

3 

•914 285 714 

13 

•477 955 737 

4 

.812 698 412 

14 

.461 474 553 

5 

.738 816 7^8 

15 

.446 5S0 277 

6 

.681 984 681 

16 

.433 055 299 

7 

.636 519 036 

17 

.420 682 290 

8 

• 599 076 740 

18 

.409 312 4 99 

9 

• 587 546 386 

19 

.398 817 307 

10 

<540 520 367 

L._.. _J 

, — 


Table 2 may be found at the end of the text. 




(d ) Aflymptotic summation of th e serie s and. amroxlmatlon curves *- To bo able to indicate 

^ V (l) > v — ( 2 ) 

approximation curves for ^7 and that pass through the computed points 

as closely as possible the limiting values for t] 0 ~ > “ mufft b© established, 

first- 

'Tn -hVil a mivTinniA I’rmnl rtRY" +.Via haiHar \ ^ 2 ^ at. n tc n f c»c annn'klnn kO\ 

" '**' Trmmmmmm jf — — “■ — — — / W^w *j '|Q ^ UWV V4>VM Tt. f 


( 2 ) 




V“ *\ / \ i- fcOv _ t / y\ 

'U** + ** ^‘2xH-l' W 


&» 


. 3L r...*? = .. . (t , 2 + f j 

0 ( n 0 2 - s 2 ) 2 % 2 - f (, c 2 - s 2 > 2 0 


74) 


^ \ 

A \ VW 

The derivative of the first scries * 0 vratten 

U- <l> /_ = 'L lkn + 3 >Wic) | (1 - ?|2p W ?) 

/ - -o - - g - r-Tl Sntl^o' 

'o 


(75) 


so that 


lim 


«n< 1 o> 


V ' 1 A in f---l (n 0 2 - l)V(lo> 


00 

H 
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since 




= \ P n (i) m 


2 + 1 
Z - 1 


A n -1 


(*> 


with 


P n- 1 « * ^rfr p n- 3 < z) + sf^rfy P n- 5 <*> + 


hence 

(z 2 - D^'iz) 


-P n (z) + (z 2 - 1 ) 


— P 1 ( z ) In - f ‘ 

2 r n \z) in z _ 1 I n . 1 


and 


1 - — 


Qn(z) 


( z) In 

~n' z - 1 


- l)Qn'(z) 


1 - (z 2 - 1 ) 


A fR^i in z.+,.l .f I 

2 ? n (z) z - 1 n-lj 


-»-l 


f n _^ is a polynomial of the (n - l) th decree in z, hence finite 


*n<D - 4 m f4i - 1) ** f4i +0 


z + 1 


Thus the series ( 75 ) can he reduced to the series ( 7*0 multiplied 
1 . ^o + 1 

O 



llm 


V * 


5 r~CO 


llm 


. V S,a ("° g + S 2) llm f-^g 

la* 1 (la 2 - £ > £=° lo* 1 V 1< 


n = n . 


'0 * 'o 

llm 


l d r—f 1 ) 


Qq that la — T" ) 

1„-1 So ^ t- 


~ t '^> — 1 

£° ^ 


(76) 


The corresponding asymptotic expression for the other series is obt aine d by differentiation 
of (74) with respect to 7] q : 




r^ 2 

w. 


+ 3 'S'g ttfl ('loK 1 ■ f) p 'an»l(^ 

. (X - t 2) hshL ^ - ("o 2 1 £ g ) H, = (1 _ {8) -&l„(l l * 3£ g ) 


fr. 2 - £ 2 } 


( n 0 2 - £ 2 ) 3 


(a) 


This derirative is also uniformly convergent, since ) converges lifce* the derivative 

n - n 0 

of a power series. Multiplication by (t] 0 2 - l) followed by another differentiation, 
gives again a uniformly convergent series 


CD 


UJ 
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L^J 


(4n + 3)(2n + l)(2n + 2)0 9lu . n (n o ) (l - C 2 )?'o_., 


[ho* + 3t S )U 0 2 - f) - 3i a (n a ♦ 3{ 2 )(1 


+ 12 f 2 n 2 fr> 2 - ii| 


n 2 _ ~ 

Ho 3 


^ '° ‘ *° 5*0 V 


As 0^^%) < g 204 ^ 2 n "+ ~ 2 rj — tlle conVer g eELCe is again readily apparent. For the 


term of the series we get 


«**< ■ (2n + 4)(2n + 3)(2n + 2) q 2 * 41 


The desired series reads then 


t^ 2 */ =5~ C“ + 3K 211 + iX 2 ^ + a )«em.i(n„) I 

— -o — 

- i d - CV^C) 
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and owing to (77) 




11m v“/ 

tj cl I dT l i “ 


Ti«Tl 


-urnim ~4y 2(1 - 5 s ) 


o ^ o *l 


IlnHlIn 8 - 3) + 3S 4 (3l„ g - 1) + gj 2 lofttyj -,■?> . lim ^ a 3*4^ 


4 - < a 


o - n 4 

TIq* 1 ° 


4 lim Q q (70) 

T] 0 =1 


hence 


i_ a_y< > / _i i* 

% ^ V*>0 „ =1 

/ e=o •-»- 


/ 


The practical interpolation curves are (fig* 7) 


a yd) / 

S4=_/ ~ 

Q o / n " n o ~ 2 ti 0 


^ln L 


1 + e3(l _T lo) 


s__r(2) 

<h/— / 3 


« 


J <*o /^o ^1 q 


H + # eo.a(i-no)j - -L [i + e o.8iT7(i-no] 


CD 

vji 


(79) 


££31 *°ll WL VCWK 
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These curves actually have the given values for tj - 1 "but 
the limiting values also coincide for Because for 


t) o ->« due to (70) ~ * and F~>1 


T1 -00 

•o 


Q 2n+1 ( T Jo^ ^ &2n+i 


TJ ®oo 

'o 


. 2n+2 


(no 2 ' 


2n+l 


(an + l)(n 0 Qg n+1 - Q^) — > (2n + 1) (a 

nn=“ \ 


- A 


2n+l ^2n+2 2n ^2n+l , 


(an + 1)A ( 


-itSJLl 


2r 


20+1 2n+2 \ 0 4n + 2 0 

& 


(an + DOgn+x n o - &L ± ^ 1 

<in+x 0 2n + 1 


o 2 

2n+l 


thus ’-2ZLL > - 

(la 2 " l )«'2n + l 1 0 =“ + 2) 


A +1 — i f (n) i — 

2' 2n 2n+2 _ 8n+3 

X % 


32“ t)q » 1 the second term in the series is already small 

compared to the first and it suffices to consider the first 
(n ■ 0): 


(*» + 3) 


Q 2 2n+1 


(tj 0 - ^Q'ai+i 


P W°> 


~ 3 — 2' 


2 L 


V* 

n®0 


2 i 0 3 (ano)‘ 


1 = - 


2 V 


Since moreover, Q 


, 0 — * — ono obtains for 5*0 
T) 0 =>« 'O 


Bn / 




1 
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The terms of the other series merely differ hy the factor 
{2n + 1} (2n + 2) , that is for n = 0 


a VC 2 ) / 

Ek_/ 

S> /’> =T io 



Translated "by J« Yanier 
National Advisory Committee 
for Aeronautics 
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TABLE 2 

COEFFICIENTS POE THE HTEEFGBOMETRIC SERIES TO COMPOTE THE SPHERICAL FUNCTIONS 


n 

“i 

°2 


“4 

°5 

a 6 

°7 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

1 6 

17 

18 

19 



0.100000000 

.071428571 

.055555556 

.045454545 

.038461538 

■033333333 

.029411765 

.026315789 

.023809524 

.021739130 

.020000000 

.01851B519 

.017241379 

.016129032 

.015151515 

.014285714 

.013513514 

.012820513 

.612195122 

0.032142857 

.017857143 

.OU363636 

.007867133 

.005769231 

.004411765 

.003482972 

.002819549 

.002329193 

.001956522 

.OOI666667 

.001436782 

.001251390 

.001099707 

.000974026 

.000868726- 

.000779626 

.000703565 

.000638117 

0.014880952 

.006764069 

.003642191 

.002185315 

.001414027 

.000967492 

.000691066 

.000510788 

.000388199 

.000301932 

.000239464 

.000193116 

.000158004 

.000130918 

.000109688 

.000092813 

.000079230 

.000068175 

.000059085 

0.008285985 
.003186917 
.001487228 
.000787356 
.000455838 
•000282 185 
.000184034 
.000125143 
.000088064 
.000063770 
.000047313 
.000035843 
.OOOO27651 
.OOOO21672 
.OOOOI7227 
.OOOOI3865 
.000011286 
.000009279 
.000007700 

O.OO51628O 6 

.001720935 

.000708620 

.000335662 

.000175823 

.000099378 

.000059627 

.000037543 

.000024597 

.000016662 

.000011613 

.000008295 

.000006053 

.000004501 

.000003403 

.000002612 

.000002031 

.000001599 

.000001273 

0.003470553 

.001020751 

.000376066 

.000161171 

.000077082 

.000040082 

.000022268 

.000013054 

•000008001 

.000005091 

.000003346 

.000002261 

.OOOOOI565 

.000001107 

.OOOOOO798 

.OOOOOO585 

.000000436 

.OOOOOO329 

0.002464384 

.000648522 

.000216174 

.000084590 

.000037220 

.000017920 

.000009269 

.000005083 

.000002927 

.000001756 

.000001092 


oo 

VO 
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as a function of t) q . 
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Figure 2.- Maximum excessive speed (u-, /U*\ on ellipsoids of the 

1 /max 

thickness ratio — = e* for various Mach numbers M = — — . 


00 
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Figure 3. - Disturbance velocities on ellipsoid of first and second 
approximation. Comparison with Lamia’s values for the sphere 
d*/t = €* = 1. 
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Figure 4. - Ratio of the excessive velocity u-^ = u - U* compressible 
to incompressible, first and second approximation. 
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Figure 5.- Cutout from figure 4: d*/t = e* = 0.2. Comparison with 

the asymptotic values. 
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Figure 6.- Ratio of the total velocities u = U* + u-^ compressible 

to incompressible as a function of the Mach number M 
(Second approximation) . 
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by thin approximation curves (calculation values +) and comparison 


with the (dashed) approximation curves. Gh\= ^ s/n = ii~. 
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Figure 8.- Maximum velocity on the ellipsoid of the thic 

for incompressible flow (M = 0). 
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